Strong Correlations Between Fluctuations and Response in Aging Transport 
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Once the problem of ensemble averaging is removed, correlations between the response of a single 
molecule to an external driving field F, with the history of fluctuations of the particle, become 
detectable. Exact analytical theory for the continuous time random walk and numerical simula- 
tions for the quenched trap model give the behaviors of the correlation between fluctuations of the 
displacement in the aging period (0, t a ), and the response to bias switched on at time t a - In par- 
ticular in the dynamical phase where the models exhibit aging we find finite correlations even in 
the asymptotic limit t a — > oo, while in the non-aging phase the correlations are zero in the same 
limit. Linear response theory gives a simple relation between these correlations and the fractional 
diffusion coefficient. 
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Originally the fluctuation-dissipation theorem was for- 
mulated by Einstein using simple Brownian motion, i.e. 
the well known Einstein relation between the diffusivity 
of an ensemble of identical particles and their mobility 
. Modern optical techniques enable experimentalists to 
track single molecules undergoing either normal diffusion, 
or non Markovian diffusion, and in some cases anomalous 
diffusion in a very wide variety of Physical and biolog- 
ical conditions 0, H 0, ■ An issue over looked so far 
are the correlations between the response of individual 
particles to an external bias and their history of fluctu- 
ations. We will soon define these correlations precisely, 
however for the time being consider chemically identi- 
cal single molecules under going a diffusion process in 
some random medium. In the time period (0,i a ), called 
the aging period, different particles may exhibit different 
fluctuations, for example if the particles sample different 
realizations of disorder in the system. Then an exter- 
nal driving field is switched on at time t a . We realize 
that the response of different particles might be corre- 
lated with with their diffusivity in the aging period. For 
example if a particle happens to sample a region in space 
with deep traps in the aging period, it is expected to 
diffusive slowly and then also respond weakly compared 
with another particle which happened to sample rela- 
tively shallow traps. However if the aging period is long 
t a — » oo and the process ergodic, we may expect that 
diffusivity of the particles in the aging period becomes 
de-correlated from their response, since the particles had 
enough time to sample the full random energy landscape 
in which they are moving. 

Here we investigate the correlations between diffusion 
(fluctuations) and mobility (response) , found when single 
particle trajectories are analyzed. We use two well known 
models:- the Montroll- Weiss continuous time random 
walk (CTRW) model 6, 7] and Bouchaud's quenched 
trap model [8j. These models are known to exhibit 
anomalous diffusion 0, Q and aging H H El El El 
i.e. the ensemble averaged response to an exter- 
nal field switched on at time t a depends on the aging time 



(see details below) . We show that in the phase where the 
models exhibit aging and even in the limit where t a — > oo 
correlations between fluctuations and response are finite, 
in complete contrast to non-aging dynamics. Besides the 
theoretical interest in correlations in aging transport, the 
existence of such correlations can be used in principle to 
predict which members of an ensemble of diffusing parti- 
cles will respond strongly (or weakly) to a driving force. 
Our work also indicates that data analysis of the aging 
fluctuation-dissipation relations must be made with care, 
since the response and the fluctuations are generally cor- 
related. 

There are many examples of systems where such corre- 
lations might become important, we mention the recent 
experiments on the diffusion of single Lad repressor pro- 
teins on DNA jjj, where a wide distribution of diffusion 
coefficients of the proteins was found. This distribution 
is likely due to the random DNA sequence the single pro- 
tein explores. Since the single molecules have widely dis- 
tributed diffusion constants, their response to an external 
field is likely correlated with their history of diffusion, e.g. 
particles with small (or large) diffusion constants have a 
weak (strong) response respectively. Other single particle 
experiments of micro-beads diffusing in actin networks 
exhibit power law waiting times and anomalous diffusion, 
very much reminiscent of the trap and CTRW models. 
Hence these single particle systems are candidates for the 
investigation of the new correlations we consider in this 
manuscript, correlations which go beyond the Einstein 
relations between the ensemble average diffusion coeffi- 
cient and mobility. 

Model 1 We consider the well known one dimensional 
CTRW on a lattice 0,0 El El- The lattice spacing is 
a and the jumps are to nearest neighbors only. Waiting 
times between jump events are independent identically 
distributed random variables with a common probability 
density function (PDF) ip(r). After waiting the particle 
has a probability 1/2 + h/2 or 1/2 — h/2 to jump to the 
right or left respectively In the aging period (0, t a ) h = 
and the particles follow an unbiased motion, while in the 
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response period t a < t the bias is < h < 1. The total 
measurement time is t = t a + t r where t r is called the 
response time. To define the response one has to define 
the field which is responsible for the bias h. For example, 
if the particle is coupled to a thermal heat bath with 
temperature T, and driven by a uniform force field F, 
standard detailed balance conditions give h = aF/2kbT, 
when h < < 1 [Tsl Il7j . We consider later the generic case 



|r(-«)| 



(i) 



when t — > oo and 0<a<l,A>0. Specific values of a 
for a wide range of physical systems and models are given 
m • F° r example for the annealed version of the trap 
model a — T/Tg < 1 9]. In this case the average waiting 
time is infinite. 

The position of the particle at time t a + t r is X = 
X a + X r , where X a (X r ) is the displacement in the 
aging (response) periods respectively. More specifically 



Xa = ESl M ' X r = YZl Ax i > Where Ax i and 

are the random jump lengths (of length a) in the 
aging (0, t a ) and response periods (t a , t a +t r ) respectively. 
While n a and n r are the random number of jumps in the 
aging and response periods respectively. 

We investigate the correlation function ((X a ) 2 X r ) 
which is a measure for the correlation between the fluc- 
tuations in the aging period {Xa) 2 , and the response to 
the driving force switched on at time t a , X r . We define 
a dimensionless fluctuation-response (FR) parameter 



FR(t a ,t r ) — 



((XgfXr) 
((Xa) 2 )(X r ) 



1, 



(2) 



which is equal zero when correlations vanish. For 
the CTRW under investigation one can show that 
{(X a ) 2 X r ) — ha 3 (n a n r ) and 



FR(t a ,t r ) 



(n a n r ) 
(n a )(n r ) 



1. 



(3) 



Thus (n a n r ), the correlations of the number of steps in 
the aging period with the number of steps in the response 
period gives a measure for the correlations between the 
fluctuations in the displacement in the aging period and 
the response to the bias. 

Let P ta- t r ( n a, Ti r ) be the probability of making n a 
jumps in the aging period and n r jumps in the re- 
sponse period. Knowledge of this function is needed 
for the calculation of the FR parameter and other high 
order correlation functions 18] which we will discussed 
in a future publication. The paths with n a (n r ) jump 
events in the aging period (response period) clearly sat- 
isfy t na <t a <t na + l and {t na +n r 

respectively, where the subscript n in t n is for the jump 
number. Hence 

Pt a ,t r (n a ,n r ) = 



< t r + t a < t na+nr + 1 )} 

(4) 

where I(x) = 1 if the event in the parenthesis is true 
otherwise it is zero. Us-usual we make use of the double 
Laplace transform t a — > u and t r — > s of Pt a ,t r (iia,n r ), 
Pu,s (n a ,n r ) = J£° dt r e' sU / °° dt a e~ uta P ta , tr . (n a , n r ) . 
For the sake of space we do not discuss the details of 
the calculation (which will be published later). We use 
the model assumption of independent identically dis- 
tributed waiting times, namely the renewal property of 
the CTRW, and find 



P u ,s (n a ,n r = 0) 



f« (u) 



i — ip{u) ii>(s) — ip{ u ) 



(•5) 



while for n r > 1 



Pu,s {j^a i Tlr) 



s (u — s) 



1 — ip (s) ip (s) — ip (u) 
(6) 

In Eqs. 1)51611 (u) and i/j(s) are Laplace transforms of 
the waiting time PDF. Note that Eqs. l r >l(il) give the 
proper normalization since J2n =o X^=o Pu,s{n a , n r ) = 
l/(us). Using Eq. JgJ and (n a n r ) u ^ s = 
^^=0 ^^=0 n an r PuA n a, n r ) wc find 



•tp (s) — ip (u) ip (u) 



s (u 



2 r 



l-ip(u) l-i/>(s) 



(7) 



and the averages (n a ) u ,s — ip(u) / {us[l — ^(u)]}, 

tp(s) - %p{u) 



s{u — s) 1 — ip(u) 1 — ip(s) 



(8) 



In principle once the double Laplace inversion of Eqs. 
(|7I8|) is made, we can calculate the FR parameter. 

If the dynamics is Markovian, namely the waiting time 
PDF is exponential ip(t) = Rexp(-Rt) 



(n a n r ) = (n a )(n r ) = Rt a Rt ri 



(9) 



and FR(£ Q , t r ) — 0. For any non-Markovian process with 
a non-exponential waiting time PDF the FR parameter 
is generally not equal zero. 

If the average waiting time (r) = J Q tiP{t)At is fi- 
nite and in the limit t a — > oo the fluctuation-response 
parameter Eq. J5J) satisfies 



lim FR(i a ,* r ) = 0, 



(10) 



and the correlations are lost in this limit. To see this use 
Eq. J7J in the u — > limit, the small u expansion %p(u) ~ 
1 - u(t) to find {XlXa) ~ ha 3 lp j^-pj^. Interestingly 
this result is valid for any s > 0, namely both for short 
and long response times t r . Hence when t a — > oo we find 
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{X 2 X r ) ~ ha 3 j^j£j and similar calculations for {(X a ) 2 ), 
and (X r ) complete the proof of Eq. H10|) . If the aging 
period is much larger than the finite time between jumps 
the response is not correlated with the fluctuations, since 
the particles had enough time to equilibrate. 

A very different behavior is found in the common sit- 
uation |j, 0, where the average waiting time is infinite 
namely when < a < 1 in Eq. Q. Then using Eq. JjJ, 
in the limit of small s and u 

1 u a - s a 



{(X a y X r ) ~ ha- 



A 2 (u - s)u 2a s 



2a „l+a ' 



(11) 



where the expansion ip(u) ~ 1 — Au a was used. Skipping 
the technical details, we analytically invert the double 
Laplace transform in Eq. to the double time domain 
and find 



J-2Q 



{(X a ) 2 X r ) ~ha 



3 "r 



A 2 ' 



(12) 



which is valid in the limit of long times t a and t r . The 
scaling function in Eq. I|12|) is a hypergeometric function 



2 F l [1, -a; 1 



„2a 



r 2 (i 



r(l + 2a) 



(13) 



Eq. (|12[1 is a main result of this paper, since it shows 
that even in the long aging time limit a non-trivial corre- 
lation between fluctuation and response exists. The hy- 
pergeometric function in Eq. I|13[) is tabulated in Math- 
ematica hence the solution is not a formal expression. 



Using ((X a y 
the model [jj 



(X r ) 



Ar(l+a) 



AT(l+a) 

fluctuation-response parameter is 



, and the aging response of 



the dimcnsionlcss 



FR(x) 



2 Fi(l,-a;l 



: r 2 (l+a) 
T(1+2q) 



(1 



1, (14) 



where x = t a /t r . If a = 1 we have FR(x) = indicating 
that the non-trivial correlations are found in the limit of 
long times, only for anomalous processes with a < 1. Eq. 
H14JI is valid in the limit t a — > oo and t r — » oo their ratio 
x remaining finite. 

Comparison between simulations of the CTRW process 
and Eq. I|14|) for a = 1/2 and a — 3/4 is made in Fig. 
^ The figure illustrates that the correlations between 
fluctuations and response becomes larger as x — t a /t r 
is increased. This is the expected behavior, the larger 
is the aging time t a compared with the response time 
t r the stronger is the correlation, since if t r >> t a the 
particle already "forgot" its behavior in the aging period. 
The Fig. also demonstrates that as a is decreased the 
correlations get stronger. In the simulation t a = 3 * 10 6 , 
■0(t) = Q f <r~( 1+Q ) for t > 1 otherwise it is zero. 

Using Eq. I|14|) we find in the limit x = t a /t r << 1 



FR(x) 



1 



r 2 (i + a) 

r(l + 2a) 



+ 0{x), (15) 




x=t IX 



FIG. 1: The fluctuation-response parameter FR(a;) as a func- 
tion of x = t a /t r . Simulations and theory agree well without 
fitting. We use a = 3/4 (crosses) and a = 1/2 (stars) and see 
that correlations are stronger when a — 1/2. 



namely weak correlations between fluctuations and re- 
sponse. In the opposite limit, of the aging regime of 
x » 1 the correlations are stronger, and we find 



FR(x) 



o|T(a)p 
T(2a) 



1 1 



1 



a x L 



(16) 



The leading term gives the non-trivial behavior of the 
fluctuation-response parameter when t a /t r — > oo. We 
find the bounds < linx^oo FR(x) < 1 where the lower 
bound with zero correlations corresponds to a — > 1 and 
the upper bound of strong correlations is found when 
a -> 0. 

Applying linear response theory to Eq. (|12|) , yields the 
connection between the correlation function and Physi- 
cally observable parameters. In this same limit aging 
Einstein relations between the ensemble average response 
and the fluctuations in the absence of the field are valid 
[Hull. We find usin g h = aF/2k b T -> 



({x a yx r 



2FD 2 a t 2 r a 
k b T ' 



(17) 



where D a — a 2 /2A is the fractional diffusion coefficient, 
which according to its definition is (X 2 ) ~ 2D a t a /F(l + 
a) 0113. Eq. 

is important since it shows that the 
transport coefficient D a and the exponent a, describ- 
ing the fluctuations in the absence of the external driving 
field, are the only system parameters needed for the de- 
termination of the correlation between fluctuations and 
the response. 

Model 2 As shown by Feigelman and Vinokur [2l| 
transport in disordered systems with quenched disorder 
may exhibit aging effects. Hence it is natural to check 



4 



- 










i 






f ..- 

* / 








■L + 




+ 



5 10 15 20 25 30 35 40 

t IX 

a r 



FIG. 2: The FR parameter for the quenched trap model. 
The + is for the high temperature phase T = 5T g /2 namely 
T > T g which does not exhibit correlations FR ~ 0, while 
the stars are for the aging phase T = T g /2 which exhibits 
non-trivial correlations. The solid curve is the CTRW theory. 

if fluctuations and response are correlated for models of 
quenched disorder, and if so how do they compare with 
those we obtained analytically in the annealed CTRW 
model? In particular do quenched models also exhibit 
a transition between an aging regime with strong cor- 
relations to a regime with vanishing correlations when 
t a — ► oo? For that aim we consider the quenched trap 
model on a one dimensional lattice 0, 0|. Each lat- 
tice site i has a fixed random energy Ej > 0, which is 
the energy barrier the particle has to cross in order to 
jump from i to i + 1 or i — 1. The energy barriers are 
all independent identically distributed random variables 
with a common PDF p(E) = {T g )- l e- E / T <> . The PDF 
of escape times from site i is exponential with a mean es- 
cape time Ti — exp(Ei/T). Notice that according to this 
Arrhenius law small fluctuations in the energy may lead 
to exponentially large fluctuations in the waiting time in 
site i. After waiting in the trap for a random time the 
particle has a probability 1/2 of jumping left or right if 
the system is not biased. It has a probability (1 ±h)/2 of 
jumping left or right when the bias is not zero. In sim- 
ulations one lets the system evolve without bias in the 
aging period, (0,t a ) and then a bias is switched on. 

It is well known that the model exhibits aging behav- 
iors when T < T g || • Bertin and Bouchaud showed that 
the aging exhibits linear and non-linear types of response, 
depending on the magnitude of h [T7J . Here we consider 
only the linear response regime of h — > 0. The trap model 
is not an exactly solvable model since the effect of the 
quenched disorder is to induce non-independent waiting 
times in the random walk. Hence we investigate the trap 
model using numerical simulations. 

Simulation result for the FR parameter are shown in 



Fig. |2 We choose a parameter set which is known to 
exhibit aging behavior 17] t a — 10 6 , h = 0.008, and 
vary t r . We choose two values of temperature T the first 
is in the high temperature ergodic phase T/T g = 5/2, 
in this case we see that the FR parameter is practically 
zero with some small deviations due to the finite time 
of simulation. For T/T g = 1/2, namely in the aging 
phase, we see strong correlations between the fluctua- 
tions and the response especially when t r < t a . We 
also plot the FR(x) parameter for the CTRW model us- 
ing the exponent a = 2T/T g /(1 + T/T g ) = 2/3 (a is 
the exponent describing the averaged response function 
[l7j)- The Fig. clearly demonstrates that the FR pa- 
rameter in the CTRW theory is smaller than the cor- 
responding FR parameter of the quenched trap model. 
In the quenched trap model, unlike the CTRW process, 
the diffusion is strongly correlated in the sense that a 
particle once returning to a specific trap will recall its 
waiting time for that trap, hence the FR parameter for 
the quenched model is larger than the one found for the 
CTRW. 

To conclude we showed that for any non-Markovian 
CTRW, correlations between the fluctuations in the ag- 
ing period and the response are finite if the aging period 
is finite. Both for the CTRW with a < 1 and for the 
quenched trap model with T < T g a non-trivial FR pa- 
rameter was found, even in the limit of large aging times. 
These correlations are found in the phase where the mod- 
els exhibit aging transport, and hence we suspect that 
generally systems which exhibit aging may exhibit simi- 
lar correlations between fluctuations and response. The 
advance of single molecule tracking makes this research 
timely, since we showed that the response of individual 
particles to an external field depends on their history of 
fluctuations. 
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